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Abstract. This article is devoted to the study of muhivalued semigroups and their as- 
ymptotic behavior, with particular attention to iterations of set- valued mappings. After 
developing a general abstract framework, we present an application to a time discretiza- 
tion of the two-dimensional Navier-Stokes equations. More precisely, we prove that the 
fully implicit Euler scheme generates a family of discrete multivalued dynamical systems, 
whose global attractors converge to the global attractor of the continuous system as the 
time-step parameter approaches zero. 

1. Introduction 

The variety of questions related to evolution equations arising from fluid mechanics prob- 
lems constitutes a challenging and fascinating area of mathematics, which has attracted 
the attention of a wide number of researchers for many years. One important aspect, 
among others, is the understanding of the behavior of solutions to differential equations 
as time goes to infinity. For autonomous systems, this translates into the study of the 
properties of a semigroup of operators {S{t)}t>o, also called a dynamical system, acting 
on a phase space X, typically a Banach space or, more generally, a complete metric space 
[H m El 123]. Notice that the parameter t could be regarded as discrete, if dealing with a 
difference equation, or as continuous, in the case of a differential equation. 

When global existence and uniqueness of solutions can be proved, dynamical systems 
arise as solution operators assigning to a certain initial condition x G X the corresponding 
solution x{t) = S{t)x of the evolution problem under consideration. Unfortunately, in 
many instances, uniqueness of solutions may be hard to prove, or even out of reach. In 
this case, one has to deal with a so-called multivalued semigroup of operators, for which 
S{t)x is the set of all possible solutions at time t. Multivalued semigroups have been 
investigated by many authors, and are particularly powerful in the study of abstract 
differential inclusions, doubly nonlinear equations, gradient flows and stochastic partial 
differential equations [21 El El [II [13 [IB] • 

From the large time behavior viewpoint, the most relevant object is the so-called global 
attractor, namely, the unique compact subset of the phase space which is at the same 
time invariant and attracting. As noted in [27j, it is of crucial importance to understand 
whether the longterm dynamics of a system possessing global attractor can be properly 
approximated by discrete attractors of discrete dynamical systems generated, for example. 
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by numerical schemes associated to the evolution problem under concern. This issue has 
been widely investigated and the interested reader is referred to, e.g., [101 [HI [121 [ISl [ISl 
[211 [26]. 



1.1. The physical model and its approximation. Let C be a bounded domain 
with smooth boundary dQ. For t > and > 0, we consider the two-dimensional 
Navier-Stokes equations [201 [22] 




dtu — uAu + {u ■ V)u + Vp = /, 
divM = 0, 



where / is an autonomous incompressible forcing term. The system is supplemented with 
the nonslip boundary condition 

(1.2) u(a;,t)U6an = 0, 

and the initial condition 



(1.3) m(x, 0) = uo(x), 

along with its time discretization provided by the fully implicit Euler scheme [23] 

(1.4) ^ - z/Am" + (m" ■ V)m" + Vp" = /, M° = Mo- 
k 

Due to its nonlinear nature, uniqueness of solution to such a numerical approximation 
can be proved only by restricting the time-step parameter A; > to be small enough. 
To be more precise, such a restriction depends on the initial datum uq and does not 
therefore allow to define a single-valued discrete semigroup of operators in the classical 
sense. Nonetheless, we will be able to show that, in fact, this difficulty can be overcome 
by defining a family of multivalued discrete semigroups {Sk, < k < ki}, where ki > 
is constant and independent of uq. 

Taking advantage of previous results contained in , we will address the issue of the 
existence of a family of discrete global attractors Ak of Sk and we will prove that Ak ^ A 
as /c — )■ in a suitable sense, where A is the global attractor of the single- valued dynamical 
system generated by fll.ip - fll.3p . 



1.2. Structure of the paper. In the next section we develop, along the lines of [H], 
the abstract machinery of multivalued dynamical systems needed to prove the results 
described above. In particular, we will focus on the existence of multivalued global at- 
tractors, and prove a convergence theorem for families of discrete attractors depending 
on a parameter. Section [3] is dedicated to the applications of the abstract theory to the 
Navier-Stokes equations and their approximation. It is shown that the fully implicit Euler 
scheme generates a multivalued discrete dynamical system whose asymptotic dynamics is 
properly related to the one of the Navier-Stokes equations. 
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2. The Abstract Framework 

We here develop in an abstract way the main tools needed for the study of multivalued 
semigroups and their asymptotic behavior. Some results are more or less already known 
thanks to the works [2] and [H] . We also refer to [B] and to the more recent work [5] for a 
fairly complete overview on the subject. Except for the approximation results, equivalent 
statements of the theorems in this section can be found in [14J. We report here alternative 
proofs for completeness and review purposes. 

While the above papers mainly dealt with continuous-time semigroup, we are more in- 
terested in discrete dynamical systems generated by numerical schemes for which unique- 
ness of solutions may not hold. In particular, we will address the question of under what 
conditions a discrete multivalued semigroup can approximate the longtime behavior of a 
continuous one. 

2.1. Multivalued Semigroups. Let {X, || • ||) be a real Banach space, and let T be 
either = [0, oo) or N. A one-parameter family of set-valued maps S{t) : 2^ — )■ 2"^ is a 
multivalued semigroup (m-semigroup) if it satisfies the following properties: 

(5.1) S{0) is the identity on 2^; 

(5.2) S{t + t) = S{t)S{T), for all t, r G T. 

To simplify the notation, if x G X we will write S(t)x in place of S{t){x} and, as 
customary, for any set i3 G 2^^, we will assume 

S{t)l3 = [j S{t)x. 

The m-semigroup is said to be closed if it fulfills the further property: 

(5.3) S(t) is a closed map for every t G T, meaning that if x„ — )■ x and y„ G S(t)xn is 
such that yn y, then y G S(t)x. 

Remark 2.1. The notion of continuity of a set- valued map is not as immediate as in the 
single-valued case. An m-semigroup is called upper semicontinuous if given x G X and a 
neighborhood U{S{t)x) of S{t)x, there exists 5 > such that 

\\x-y\\<5 S{t)y cU{S{t)x). 

On the other hand, S{t) is defined to be lower semicontinuous if given Xn ^ x and 
y G S(t)x, there exists yn G S(t)xn such that yn y- Finally, S(t) is continuous if it 
is at the same time lower and upper semicontinuous. In general, it is not true that a 
continuous m-semigroup is closed (as it is in the single- valued case), the problem being 
that S(t)x might not be a closed set (if S(t)x is closed for any x E X, the S{t) is said to 
have closed values). Nonetheless, if S(t) is upper semicontinuous and has closed values, 
then S(t) is closed (see pQ). 

The positive orbit of B, starting at t G T, is the set 

T>t 

where we agree to set 7(i3) = ■joiB). A function ?/ : T — )■ X is said to be a trajectory 
starting at ?/o G X if y{0) = yo and y{t + r) G S{t)y{T) for every t, r G T. 
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Limit Sets. For any B E 2^ , the set 

is called the u-limit set of B. The following characterization of u{B) holds true in the 
multivalued case, and the proof is identical to the single- valued counterpart. 

Lemma 2.2. A point x E X belongs to uj{B) if and only if there are sequences t„ — oo 
and Xn G S{tn)B such that x„ — ?■ x as n — t- oo. 

A nonempty set G 2"^ is invariant for S{t) if 

S{t)B = B, Vt G T. 
If S{t)B C B, then B is said to be positively invariant. 

Dissipativity. A set Bq G 2^ is an absorbing set for the m-semigroup S{t) if for every 
bounded set i3 G 2"^ there exists te G T such that 

S{t)B C Bo, Vt > tB. 

Given two nonempty sets B,C G 2^ , the Hausdorff semidistance between B and C (in X) 
is defined as 

dist{B,C) = supinf ||6 — c||. 

Notice that, in general, dist(i3,C) 7^ dist(C,i3). A nonempty set C G 2"^ is attracting if 
for every bounded set B we have 

lim dist{S {t)B,C) = 0. 

t—^oo 

The m-semigroup S{t) is called dissipative if it possesses a bounded absorbing set. Some 
authors prefer to require compactness (instead of only boundedness) in the notion of 
dissipativity, a strategy which has been successful expecially in parabolic problems [16] . 



The Global Attractor. A nonempty compact set ^ G 2"^ is said to be the global 
attractor of Sit) if 

(A.l) ^ is invariant; 

(A. 2) A is an attracting set. 

Remark 2.3. The global attractor, if it exists, is necessarily unique. Moreover, it enjoys 
the following maximality and minimality properties: 

(i) let ^ be a bounded set satisfying (A.l) Then Ad A; 

(ii) let ^ be a closed set satisfying (A. 2) Then A C A. 

In order to state the result on the existence of the global attractor, we need a definition. 
Given a bounded set i3 G 2^, the Kuratowski measure of noncompacteness a{B) of B is 
defined as 

a{B) = inf {5 : B has a finite cover by balls of X of diameter less than 5}. 

We list hereafter some properties of a. 
(K.l) a{B) = a(B); 
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(K.2) Bi C B2 implies that a{Bi) < ^(-Ba); 
(K.3) a{B) = if and only if B is compact; 

(K.4) if {Bt}t£T is a family of nonempty closed sets such that Bt^ C Bt^ for t2 > ti and 
lim4_j.oo a('Bt) = 0, then B = HteT nonempty and compact; 

(K.5) if {Bt}teT and B are as above, given any t„ — )■ cx3 and any Xn € Bt„, there exist 
X & B and a subsequence x^^. — )■ x. 

Theorem 2.4. Suppose that the dosed m-semigroup S{t) possesses a bounded absorbing 
set Bo e 2^ and 

(2.1) lim a{S{t)Bo) = 0. 

t-+oo 

Then u){Bo) is the global attractor of S{t). 

Proof. First, we prove that oj^Bq) is nonempty and compact. Let to ^ T be such that 

S{r)Bo C Bo, Vr > to- 

If t > to, we have 

T>t tGT 

= U 5(t - to)5(r + to)Bo clJSit- to)Bo = S{t - to)Bo. 



Therefore, in hght of (K.l) (K.2) and (EH), 



lim Q;(7t(i3o)) = lim a(7t(i3o)) = 0. 
Since the sets 7t(-Bo) are nested, from (K.4)| we conclude that 

uj{Bo) = {yiM) 

teT 

is nonempty and compact. 

To prove uj{Bo) is attracting, argue by contradiction, and assume that there exist e > 
and sequences t„ — 00 and ?/„ G S{tn)Bo such that 

inf \\yn - y\\ > e. 



Since yn G S{tn)Bo, it follows that G 7t„(i3o). By properties (K.4) and (K.5) , there exist 
z G uj{Bo) and a subsequence yn^. G S'(tnj,)So such that yn^ z, which is a contradiction. 

It remains to show that uj{Bo) is invariant. Let t„ — )■ 00 and ?/„ G 5'(t„);Bo be any 
sequences. We claim that there exists y G uj{Bo) such that ?/„ — )• y up to a subsequence. 
Since ^(^Bo) is attracting, we know that 

lim dist(^(t„)i3o,w(i3o)) = 0. 



As a consequence. 



lim inf — -^H = 0, 



from which we deduce the existence of a sequence Zn G uj{Bo) such that 

lim - Zn\\ = 0. 
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From the compactness of uj{Bo) we readily get y G uj{Bo) and a subsequence n^. such that 
y, and, in turn, y^^ y. 
Let X G oj^Bo) and t G T. By Lemma 12.21 there exist sequences t„ — )■ oo and G 
S(tn)Bo such that ?/„—)■ x as n oo. Since 

5(t„)i3o = 5(t)S(t„-t)i3o, 

we obtain that yn G S{t)zn, with 2^ G S'(t„ — t)Bo. By the above claim, there exists 
y G w(i3o) such that z„ — )■ y up to a subsequence. The fact that S{t) is closed then yields 
the following implication 

Zn^y, S{t)zn3yn^x xeS{t)y, 

namely the inclusion uj{Bo) C S{t)uj{BQ). Turning to the opposite one, we now know that 

S{t)co{Bo) C S{r)S{t)uj{Bo) = S{t + r)uj{Bo) 

for every t, r G T. Since u{Bo) is attracting, for every neighborhood U{u{Bo)) of u{Bo), 
there exists G T such that 

S{t + t)uj{Bo) C W(a;(So)), Vr > t^,. 

Thus S{t)uj{Bo) C W(c(;(i3o))) where U{u{Bo)) is arbitrary, so that 



S{t)u{Bo) C w(i3o) = u;(i3o), Vt G T, 
which concludes the proof. □ 

2.2. Discrete Approximation Of Multivalued Semigroups. We now turn our at- 
tention to the approximation of continuous-time m-semigroups by discrete ones, exploring 
how certain properties that hold true in the discrete regime are carried over to the con- 
tinuous limit. 

Remark 2.5. Given a set- valued map S* : 2"^ — 2"^, we can define a discrete m-semigroup 
by 

S{n) = S", yn G N. 



Properties (S.l) and (S.2) are trivially satisfied, and, in this case, we will write that 
{'S'jneN (instead of the redundant {S""'}„gN) is a discrete m-semigroup. 

Turning to our approximation problem, let {S'(i(:)}jgiR+ be a closed m-semigroup, Kq a 
positive constant, and consider a family of discrete closed m-semigroups {5"^, < k < 
^ojnGN where, for each fixed k, the map Sk : 2^ ^ 2^ satisfies the usual semigroup 
properties 

Si = ld,x , Sr'" = S^S^, Vn, m G N. 

Such maps arise in the study of numerical schemes associated to evolutionary equations, 
in which either the discretization problem or the differential system (or both) might not 
enjoy any uniqueness property of solutions. From the point of view of the longtime 
behavior of solutions, it is therefore interesting to understand under what conditions the 
asymptotic features of the continuous-time m-semigroup can be properly approximated 
by the discrete ones. 
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Given two nonempty sets i3, C G 2"^, we write 



B-C = {b-c: beB,ceC} 



and 



= sup 



The following theorem is a generalization to m-semigroups of a result proven in 

Theorem 2.6. Let S{t) be a closed m- semigroup, possessing the global attractor A, and 
let {Sk, < k < KojneN be a family of discrete closed m-semigroups, with global attractor 
Ak- Assume the following: 

(HI) there exists ki G (0, kq] such that the set 



(2.2) 



IC= [j Ak 



IS bounded in X ; 
(H2) there exists to > such that for any T* > to, 



(2.3) 



lim sup \\SJ^x — S{nk)x\ 

''^^ x£Ak,nke[to,T*] 



0. 



Then, 
(2.4) 



limdist(A,^) = 0. 



Proof. Let e > and k G (0,ki]. Since K, is bounded and A is attracting, there exists 
te > to > such that 



dist{S{t)IC,A) < \/t>t,. 



Let now Xk & Ak, and pick 



nk 



te + l 
k 



Thanks to the invariance of Ak under Sk, there exists yk G Ak such that Xk G S^'^yk-, and 



by (H2) 



\\xk - S{nkk)yk\\ < \\S^''yk - S{nkk)yk\\ < -, Wk < k^, 
with a proper choice of > 0. As a consequence, if y G S{nkk)yk, 

dist(^fc,^)= sup inf||xfc— x||< sup — ?/|| + dist(?/, ^) 

< sup \xk- S{nkk)yk\\ ^\^t{S{nkk)yk,A) 

< sup \xk — S(nkk^yk\ + dist(5'(nfc/i;)/C, ^) < e, 

and the theorem is proved. 



□ 
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3. Applications to the Navier-Stokes Equations 

In this section, we apply the above abstract framework to the fully implicit Euler approx- 
imation of the Navier-Stokes equations. It is known that such a scheme does not generate 
a single- valued discrete semigroup, since uniqueness of solutions holds under a restriction 
on the time-step parameter which depends on the initial datum. To circumvent this diffi- 
culty, we will show that such a scheme generates a family of closed discrete m-semigroups 
depending on the time-step parameter, whose related attractors converge to the attractor 
of the dynamical system generated by the Navier-Stokes equations. 

3.1. Function spaces. Let C be a bounded domain with smooth boundary dQ. 
For p e [l,oo] and A; e N, we denote by LP(fi) = {LP{Q)y, W^iQ) = {H''{Q)Y, and 
Hg(i7) = {Hq{Q)Y the usual Lebesgue and Sobolev spaces of vector- valued functions on 
Q. Setting 

V= {ue C^{n,M.'^) : divM = 0}, 
we consider the usual Hilbert spaces associated with the Navier-Stokes equations 

H = closure of V in L^(f2), 
V = closure of P in 

where we denote by | ■ |, (■, ■) and || ■ ||, ((■,■)) the norm and the scalar product in H and 
in V, respectively. Also, we indicate by V* the dual space of V, endowed with the usual 
dual norm || ■ ||*, and by (■, ■) the duality pairing between V and V*. Calling 

P : = H®H^ H 

the Leray orthogonal projection, the Stokes operator is defined as 

A = -PA, dom(A) = u^in) n V. 

It is well known that the operator A is self-adjoint and strictly positive. Moreover, 
dom(y4^/2) = V and 

\\u\\ = \Vu\ = \A^^\\, \/ueV. 

Setting 

B{u,v) = P[{u-V)v], 
system fll.ll) - fll.3p can be rewritten as an abstract evolution equation of the form 

^2 (u + uAu + B{u,u) = f, 

\m(0) = uo, 

where / = Pf, since we are assuming incompressible forcing. As proved in [23], problem 
(13. ip generates a continuous and dissipative single-valued dynamical system S(t) : H — )■ 
H, which possesses the global attractor A, bounded in V. 
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3.2. The implicit Euler scheme. A possible time discretization of (13. ip is given by the 
fully implicit Euler scheme 

(3.2) + z/Au" + m") = /, M° = uo, 

k 

where n > 1 and /c > is the time-step. Our goal in this section is to prove that (13. 2p 
generates a closed discrete m-semigroup {S'jtjngN. Indeed, as in the stationary Navier- 
Stokes problem, a solution to (13. 2p is not in general unique. In particular, uniqueness of 
solutions may be proved requiring k to be bounded from above by a constant depending 
on the initial datum. More specifically, such a bound vanishes as the norm of uq tends to 
infinity, and this makes it impossible to properly define a single-valued semigroup acting 
on the whole phase space. Nonetheless, by means of the tools devised in Section [21 we will 
still be able to give a description of the longtime behavior of the discretized Navier-Stokes 
equations, and discuss its convergence to the time-continuous asymptotic dynamics. 

Notation. Throughout the section, C and Q(-) will denote a generic positive constant 
and a generic increasing positive function, respectively, whose value may change even 
in the same line of a certain equation. Unless otherwise stated, these quantities will be 
independent of fc,n and of the initial datum uq. In general, they might depend on the 
structural quantities of the system (z/, /, VL). 

Fix A; > 0. For w E H, we look at the problem 

(3.3) u + kuAu + kB{u, u) = w + kf, 
for which we seek solutions in the following weak sense. 
Definition 3.1. A vector u E V is a solution to (13.31) if 

(3.4) {u,v) + kh'{{u,v)) + kb{u,u,v) = {w,v) + k{f,v), Wv E V, 

where b{u, v, w) = {B{u, v),w) is the usual trilinear form associated to the weak formula- 
tion of the Navier-Stokes equations. 

Remark 3.2. Recall that the trilinear form b satisfies the following properties: 

(3.5) \b{u,v,w)\<C\u\^/^\Au\^/^\\v\\\w\, W E dom{A),v E V,w E H, 

(3.6) \b{u,v,w)\<C\u\^/^\\u\\^/^\\v\\\w\^/^\\w\\^/\ W,v,weV, 

(3.7) b{u,v,v) = 0, \fu,vEV. 

It is a classical result that a (possibly not unique) solution to (13. 3p exists. Moreover, 
any solution u E V satisfies the energy estimate 

(3.8) \u\^ + ku\\u\\^ <\w\^ + Ck\f\\ 

where C > does not depend on k. For every w E H, define the multivalued map 
Sk:2^ ^ 2^ by 

SkW = {u E V : u solves (13. 3p with time-step k}. 

Notice that, in light of (13. Sp . the set SkW is bounded in V and therefore relatively compact 
in H, thanks to the compactness of the embedding V ^ H. 
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3.3. The discrete m-semigroup. Let us consider the discrete m-semigroup {S'^jngN 
generated by 5*^. It is now clear that, for every n G N, 

= {m" G : solves Q with time-step k}. 

Also, a vector G SJ^Uq if and only if there exists a sequence of elements u^, . . . , u^' 
such that G S^u^^^ for every z = 1, . . . ,n and u° = Uq. From (13. Sp . we infer that any 
G SkU^^^ satisfies the energy estimate 

\u'\^ + ku\\u%^ < \u'-^\^ + Ck\f\^, 

and, inductively, any G 5^?^'' fulfills the bound 

n 

(3.9) |M"|2 + A;i/5^||M*f < \uy + Ckn\f\\ 

We then have the following theorem. 

Theorem 3.3. The multivalued map associated to the implicit Euler scheme (13. 2p 
generates a closed discrete m-semigroup {Sk}n^^. 



Proof. Since properties (S.l) (8.2) are satisfied by definition, all we need to prove is that 
S*^ is a closed multivalued map for each ra G N. As j — ?■ oo, let u'^ m H and 

G 5'^M° with -)■ in H. We have to show that m" G 5'^'U°. 

Since G there exists a sequence (m°,m], . . . where G Sku'^'^ is a 

solution to 

(3.10) («}, t;) + kv{{u% v)) + A;6(m}, m}, v) = k{f, v) + {uf\ v), Wv e V. 

Also, the fact that 'm H implies the existence of a positive number M such that 

sup|M°p < M. 
i 

In view of (13. 9p . we obtain the bound 



£=1 

Thus, for every i = l,...,n, we have the following convergences (up to not relabeled 
subsequences) as j — )■ cxd: 

— )■ M*, strongly in H and weakly in V. 

Now, passing to the limit in (I3.10p . we readily get that 

{u\v) + kiy{{u\v)) + kb{u\u\v) = k{f , v) + {u'-\ v) , Vf G V. 

As a consequence, G SkU^~^ for each z = 1, . . . ,n. But then, m" G SkU^~^ C so 
5^ is a closed map for every n G M. □ 

From the energy estimate (13. 9p and the fact that a closed map has necessarily closed 
values, we have the following straightforward consequence. 

Corollary 3.4. The discrete m-semigroup {Sk}n£n has compact values, namely, the set 
S^Uq is compact in H for every uq E H. 
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3.4. Earlier Contributions. To continue our study, we first collect some results ob- 
tained in [22]. The first two concern initial data Uq G H. There exists hq > 0, independent 
of Mo, n, fc, such that the following hold: 

(D.l) For every k > 0, 

|5>o|' < (1 + fc)"lMoP + Vn > 0. 



(D.2) Let k G (0,fi;o]. There exists a constant Rq > with the following property: for 
every R> 0, there exists to = to{R) > such that 

|5'fcMo| < -Ro, Vnfc > to, 

whenever |mo| < R- Both Rq and to can be explicitly computed and do not depend 
on n and k. In other words, the set 

Bo = {v eH : \v\< Rq} 

is a bounded absorbing set for {5'A;}„eN- 

We now turn our attention to initial data uq E V and recall the main result derived in 
[25] . tailored to our case. 

Theorem 3.5. Suppose \\uo\\ < R, and let the time-step k he such that 

(3.11) k < K,{R) = min |/to, c]. 
Then the estimate 

(3.12) WSluoW < Q{R) 
holds true for every n > 1. 

Remark 3.6. Since the issue of non-uniqueness of solutions to ( 13. 3p is the main motivation 
for this work, let us briefly discuss one way to recover uniqueness of solution to f l3.2p . in the 
sense of (13. 4p with u = m" and w = u^~^. Let and u'2 be two solutions corresponding 
to the same initial data Uq G let i? > be such that ||no|| < R, and set = ~ u^. 
If k < Kir{R), from (13. 4p and (13.70 we learn that 

\u'f + fcz/||M"f = -A;6(u", u"). 

Now, using (13.61) and Theorem 13. 5[ for any n > 1 we obtain 

+ A;i/||M"f < CA;|m"|||m"||||m^|| < Q(i?)A;|M"|||M"|| 



Therefore, if we require 



we can conclude that 



k < min |fi;^(i?). 



2Q{Ry 
+ A;z/||m"|P < 0. 
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This estimate clearly yields uniqueness of solutions. However, two main drawbacks arise. 
Firstly, we necessarily need Uq & V, which rules out the possibility of defining a single- 
valued semigroup on the natural phase space H of weak solutions. Secondly, the restriction 
on k depends on ||mo||, and thus uniqueness of solutions depends, in the end, on the single 
trajectory chosen and not uniformly with respect to the initial datum. This is why we 
decided to tackle this problem exploiting the machinery of multivalued semigroups. 

3.5. The discrete global attractors. The higher order estimate in fl3.12p is not enough 
to conclude the existence of the global attractor for the discrete m-semigroup S^, since it 
requires the initial data Uq to be in the more regular space V. Moreover, fl3.1ip shows a 
dependence of on the initial data, which turns out to be unsatisfactory in the approx- 
imation process devised in Theorem 12.61 We now show how to overcome this difficulty. 
First of all, thanks to the existence of a bounded absorbing set, it is natural to consider 
only initial data uq & Bo- 
Lemma 3.7. Let k G (0, Kq], and consider a sequence u^, . . .), where m° = G Bq 
and G Slu^. Then, there exists G N such that 

(3.13) \\u^''\\<R^, 
where R^, > does not depend on n, k and Uq. 
Proof. Let 



Estimate (13.91) immediately implies 



+ 1. 



kuJ2h'f < Rl + Ckuklfl"^. 



1=1 

Arguing by contradiction, we infer that there exists G {1, . . . , n^} such that 

kukuWu^'^f < Rl + Cknk\f\^. 

Hence, 

||«4f<_^ + C'|/p. 
Since kuk > 1, the proof ends by setting 

Rl = ^ + C\f\\ 



□ 



Having in mind condition fl3.1ip . we now fix Hi = k,^,{R^) < kq. Combining together 
fl3.12p and the above Lemma 13.71 we obtain the following. 

Corollary 3.8. Let k G (0, ni], and consider a sequence (m°,m^, . . .), where = uq E Bq 
and G S\u^ . Then, there exists G N such that 

(3.14) < Vn > 1. 
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As a consequence, for every k G (0, ki], there exists eN such that 
(3.15) ||5^'=+"mo|| < Vn>l. 

Proof. The constant Ri = Q{R*) is given by fl3.12p . and it is clearly independent of n, k 
and Uq. In the first estimate in the statement of the corollary, might depend on the 
particular sequence (m°, u^, . . .) which originates from = uq. As in the proof of Lemma 
13. 7[ the choice 

1 

.k. 



+ 1, 



together with fl3.12p and the fact that < Uk, takes away this dependence and allows to 
obtain the uniform estimate fl3.15p . □ 

Notice that fl3.15p can be rewritten in the equivalent way 

ll'S'^uoll < -Ri, Vra > rifc + 1, 

and from the definition of n^, we have that, in particular, 

(3.16) ||^>o||<^i, \fnk>l + 2Ki. 

We summarize the above observations in the next theorem, which, in fact, improves the 
results in [25] . 

Theorem 3.9. Let ki > as above and k G (0, There exists a constant Ri > with 
the following property: for every R> 0, there exists ti = ti{R) > such that 

(3.17) \\SkUo\\<Ri, \fnk>ti, 

whenever \uo\ < R. Both Ri and ti can be explicitly computed and do not depend on n 
and k. Hence, the set 

(3.18) Bi = {v eV : \\v\\ < Ri} 
is a V -bounded absorbing set for {Sk}nm- 



Proof. Fix k G (0, ki\ and let \uq\ < R. From (D.2), there exists to = to{R) such that 

SIuq G So, 'ink > to. 
In view of fl3.16l) . setting ti = ti(-R) = to(-R) + 1 + 2ki, we obtain 

5>o G Bi, \/nk > ti, 
concluding the proof. □ 

Remark 3.10. The result contained in the above Theorem 13.91 improves the one in [25] 
in two directions. On one hand, we only require the initial data to be in H. This shows 
a regularization property analogous to the one enjoyed by the solution to the Navier- 
Stokes equations. On the other hand, we obtain a uniform restriction on the time-step k, 
independent of the initial data. 

Thanks to the above results, the discrete m-semigroup {5'fc}„eN satisfies the assumptions 
of Theorem \2M 

Proposition 3.11. For every k G (0,ki], there exists the global attractor Ak of the m- 
semigroup {Sk}n&N- 
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Remark 3.12. The global attractor Ak being the smallest closed attracting set of the 
phase space, from f l3.17p we obtain the inclusion 

Ak c Bi. 

The fact that Ri does not depend on k then yields 

(3.19) U AkCBi. 

Hence, the attractors Ak enjoy a uniform regularity property. 

3.6. The attractor approximation. In this paragraph, we prove that the longterm 
behavior of the semigroup S{t) generated by the Navier-Stokes equations (13. ip is approx- 
imated, in the sense of Theorem 12. 6^ by that of the discrete m-semigroup related to the 
fully implicit Euler scheme fl3.2p . The main result of this section reads as follows. 

Theorem 3.13. The family of attractors {Ak}ke{o,Ki] converges, as k ^ 0, to A, namely, 

(3.20) limdist(A,^) = 0, 

k^O 

where dist denotes the Hausdorff semidistance in H. 

Our goal is to apply Theorem l2.6l to obtain the convergence of the discrete attractors Ak 
to the continuous time attractor A of the semigroup S{t). By virtue of f l3.19p . assumption 



(HI) is automatically satisfied. The remaining of the section is devoted to the verification 



of the uniform convergence required by (H2) Since (H2) involves an estimate in terms of 



initial data belonging to Ak, till the end of the section we will assume 

Mo e Bi. 

Define K2 = minj/ti, K^(i?i)}, where is given by (13. lip . Notice that since Bi is absorb- 
ing, in view of Theorem 13.51 we have the following uniform estimate 

(3.21) sup sup||^>o|| < C. 

fcG(0A2] n>0 

which, in turn, implies 

rn 

(3.22) ^||u"-u"-i2<C'fc(m-i + l) + C, Vi = 1, 



where G Sku"^ ^ for every n = i, . . . ,m. Indeed, multiplying (13. 2p by 2kAu", we obtain 
Wu'^f - + - u'^^f + 2vk\Au^\^ 

+ 2fc6(M", m", AvJ") = 2k{f, Au""). 

Estimating the trilinear form using (13.50 and the bound (I3.2ip . we infer that 
2A;6(m", u", Am") < 2Ck\u''\^/^\\u''\\\Au"f'^ < y|v4M"p + CA;. 

Also, 

2fc(/,v4M") < 2fc|/PM"| < ^|AM"p + Cfc. 
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Thus, 

\\u"'f - + - u^'-^f + vk\AvJ'\^ < Ck. 

Now, summing over n = i, . . . ,m and neglecting the positive term IAm"^, we get 



m 



\un\^-\\u' 



Hence, a further apphcation of (I3.2ip entails (13.221) . 

For any > 0, we define the piecewise constant and the piecewise linear functions 

Uk{t) = m", t e [{n - l)k, nk) 

and 

Ufc(t) = + - te\{n- l)k, nk). 

K 

Notice that Uk{nk) = m". Also, it is easily seen that Uk solves 

(3.23) Uk + uAuk + B{uk,Uk) = f + '^k, 
where 

*fc(t) = uAiukit) - Ukit)) + B{ukit), Uk{t)) - B{uk{t),Uk{t)). 
Lemma 3.14. For any T* > 0, G ^^(CT^ V*) and 

(3.24) ll*fc|li^(o,T*;y-) < mT"). 

Proof. Let v & V he such that \\v\\ < 1, and let t G [{n — l)k,nk) be fixed. In light of 
fl3.6p . we have 

\{B{uk,Uk) - B{uk,Uk),v)\ = \b{uk,Uk - Uk,v) + b{uk - Uk,Uk,v)\ 

< C{\\uk\\ + llMfelDII-Sfc - Uk\\. 

Since the uniform bound f l3.2ip implies 

\\uk{t)\\ = <c 

and 



t — nk 



k 



k"|| + ik""M|)<a 



we infer that 

\{B{uk,Uk) - B{uk,Uk),v)\ < C\\uk - UkW < C||m" - 

Clearly, 

\{A{uk-Uk),v)\ < \\uk-Uk\\ < C||m"-m''-^||. 
Hence, for t G [(n — l)k,nk), we can conclude that 

\\-^k{t)\U<C\\u"-u"-'\\. 
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Thus, setting N* = [T*/k\, by the above bound and fl3.22p we finally get 
ll^^lli^(o,T*;y*) = / II^^WII*dt< Yl / ll^^WII'dt 

^0 J(n-l)k 
N* + l 

<CkJ2 \K - ^""i^ < kQ{T*). 

n=l 

Thus, the lemma is proved. □ 



We are now ready to verify (a slightly stronger version of) assumption (H2) for our 
discrete m-semigroup, which will conclude the proof of Theorem 13.131 

Lemma 3.15. For any T* > 0, 

(3.25) lim sup I'S'^^^o ~ S{nk)uo\ = 0. 

''^^uo€Bi,nke[0,T*] 

Proof. Let u = u{t) = S(t)uo be the solution to (13.11) . As shown in [23], S(t) satisfies the 
uniform energy estimate 

(3.26) sup sup \\S{t)uo\\ < C. 

For Uk defined as above and k G (0, k,2], consider the difference Vk = u — Uk, which is a 
solution to 

Vk + fAvk + B{vk, u) + B{uk, Vk) = 
Testing the above equation by Vk, one obtains 

1 d_ 

2dt' 

By (13. 6 p and (13. 26 p . the trilinear form can be estimated as 



l-yfcP + i^WvkW^ + b{vk,u,Vk) = -(^fc, Wfc). 



\b{vk,u,Vk)\ < C\vk\\\vk\\\\u\\ < -\\vk\\^ + C\vk\ 



and from the obvious bound 



■■2 , r^lhT, l|2 



-{^k,vk) < -ii^^feii 

we derive the differential inequality 

^Avk\'<C\vk\' + C\\^k\\l 
at 

Since Vk{0) = 0, an application of the Gronwall inequality together with Lemma 13.141 
gives 

\vkit)\' < Ce^* r ||vl>fc(s)||^ds < Ce^^l^fc||i2(o,T*;y.) < kQiT^. 
Jo 

As a consequence, 

\u{t)-Uk{t)\'<kQiT*). 
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But then 

lim sup \SkUQ — S{nk)uo\ 

'^'^^ uoeBi,nfee[0,T*] 

= lim sup sup \u" — u{nk)\ 

— lim sup sup \uk{nk) — u{nk) \ = 0, 

'^'^'^ uo eBi , nfee [0,T*] eS^uo 

and the proof is over. □ 

Acknowledgments. This work was partially supported by the National Science Founda- 
tion under the grant NSF-DMS-0906440, and by the Research Fund of Indiana University. 



References 

[1] J-P. Aubin, H. Frankowska, Set Valued Analysis, Birkhauscr, Boston, 1990. 

[2] J.M. Ball, Continuity properties and global attractors of generalized semiflows and the Navier-Stokes 
equations, 3. Nonlinear Sci. 7 (1997), 475-502. 

[3] V. Barbu, Nonlinear semigroups and differential equations in Banach spaces, Noordhoff International 
Publishing, Leiden, 1976. 

[4] A.V. Babin, M.I. Vishik, Attractors of evolution equations, North-Holland, Amsterdam, 1992. 

[5] T. Caraballo, J. A. Langa, V.S. Melnik, J. Valero, Fullback attractors of nonautonomous and sto- 
chastic multivalued dynamical systems, Set- Valued Anal. 2 (2003), 153-201. 

[6] T. Caraballo, P. Marm-Rubio, J.C. Robinson, A comparison between two theories for multi-valued 
semiflows and their asymptotic behaviour. Set- Valued Anal. 11 (2003), 297-322. 

[7] V.V. Chepyzhov, M.I. Vishik, Attractors for equations of mathematical physics, Amer. Math. Soc, 
Providence, 2002. 

[8] C.B. Gentile, J. Simsen, On attractors for multivalued semigroups defined by generalized semiflows, 

Set- Valued Anal. 16 (2008), 105-124. 
[9] J.K. Hale, Asymptotic behavior of dissipative systems, Amer. Math. Soc, Providence, 1988. 

[10] J.K. Hale, X. Lin, G. Raugel, Upper semicontinuity of attractors for approximations of semigroups 
and partial differential equations. Math. Comp. 50 (1988), 89-123. 

[11] E. Hill, A.T. Siili, Approximation of the global attractor for the incompressible Navier-Stokes equa- 
tions, IMA J. Numer. Anal. 20 (2000), 633 667. 

[12] N. Ju, On the global stability of a temporal discretization scheme for the Navier-Stokes equations, 
IMA J. Numer. Anal. 22 (2002), 577-597. 

[13] A.V. Kapustian, J. Valero, Attractors of multivalued semiflows generated by differential inclusions 
and their approximations, Abstr. Appl. Anal. 5 (2000), 33-46. 

[14] V.S. Melnik, J. Valero, On attractors of multivalued semi-flows and differential inclusion. Set- Valued 
Anal. 6 (1998), 83-111. 

[15] V. Pata, S. Zelik, A result on the existence of global attractors for semigroups of closed operators, 

Commun. Pure Appl. Anal. 6 (2007), 481-486. 
[16] J.C. Robinson, Infinite- dimensional dynamical systems, Cambridge University Press, Cambridge, 

2001. 

[17] R. Rossi, S. Segatti, U. Stefanelli, Attractors for gradient flows of nonconvex functional and appli- 
cations. Arch. Ration. Mcch. Anal. 187 (2008), 91 135. 

[18] G. Schimperna, S. Segatti, U. Stefanelli, Well-posedness and long-time behavior for a class of doubly 
nonlinear equations. Discrete Contin. Dyn. Syst. 18 (2007), 15-38. 

[19] J. Shen, Convergence of approximate attractors for a fully discrete system for reaction- diffusion 
equations, Numer. Funct. Anal, and Optim. 10 (1989), 1213-1234. 

[20] H. Sohr, The Navier-Stokes equations, Birkhauser, Basel, 2001. 



18 



M. COTI ZELATI, F. TONE 



[21] A.M. Stuart and A.R. Humphries, Dynamical systems and numerical analysis, Cambridge University 

Press, Cambridge, 1996. 

[22] R. Temam, Navier-Stokes equations and nonlinear functional analysis, CBMS-NSF Regional confer- 
ence Series in Applied Mathematics, SIAM, Philadelphia, 1983. 

[23] R. Temam, Infinite- dimensional dynamical systems in mechanics and physics, Springer-Verlag, New 
York, 1997. 

[24] R. Temam, Navier-Stokes equations, AMS Chelsea Publishing, Providence, 2001. 

[25] F. Tone, D. Wirosoctisno, On the long-time stability of the implicit Euler scheme for the two- 
dimensional Navier-Stokes equations, SIAM J. Numer. Anal. 44 (2006), 29-40. 

[26] F. Tone, X. Wang, Approximation of the stationary statistical properties of the dynamical system 
generated by the two-dimensional Rayleigh-Benard convection problem. Anal. Appl. 9 (2011), 421- 
446. 

[27] X. Wang, Approximation of stationary statistical properties of dissipative dynamical systems: time 
discretization. Math. Comp. 79 (2010), 259-280. 

Indiana University - Department of Mathematics 
Rawles Hall, Bloomington, IN 47405, USA 
E-mail address: micotize@indiana.edu (M. Coti Zelati) 

University of West Florida - Department of Mathematics and Statistics 

Pensacola, FL 32514, USA 

E-mail address: ftone@uwf.edu (F. Tone) 



